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This paper gives a brief review on the recent work on fractional Klein-Gordon 
fields, in particular on the Casimir effect associated to fractional Klein-Gordon 
fields in various geometries and boundary conditions. New results on Casimir 
piston due to a fractional Klein-Gordon massive field arc given. 



1. Introduction 

The possibility of the fractal structure of spacetime has gained ground lately 
due to the work on quantum theory of gravity. One promising theory known 
as causal dynamical triangulations proposed by Ambj0rn, Jurkiewicz and 
Loll 1-3 is a non-perturbative approach to quantum gravity which shows 
how the spacetime has fractal structure, with the number of dimensions 
changes smoothly from the usual four at large scales to two at the Planck 
scale. Another approach of quantum gravity, the asymptotically safe quan- 
tum Einstein gravity 4,5 also has a similar property such that spacetime 
evolves from a lower dimensional scale-invariant structure at high energy 
or micro scales to that of normal four-dimensional geometry at large scales 
or low energy. In these theories, spacetime geometry cannot be understood 
in terms of a single metric. There is a need to introduce a different effec- 
tive metric at each energy scale. Modesto 6 has applied a similar technique 
used in asymptotically safe quantum gravity to loop quantum gravity and 
obtained a similar fractal-like spacetime structure. Numerical simulations 
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of the propagation of a scalar particle in a dynamically triangulated space- 
time, with a discretized Einstein dynamics have been carried out. It is found 
that the spectral dimension of the microscopic spacetime is two and it be- 
comes four for large scales 4,7-10 . Another proposal of a Lorentz invariant 
theory in a fractal spacetime requires the spectral dimension and the ul- 
traviolet Hausdorff dimension of the spacetime to be both equal to two 11 . 
The fractal nature of quantum spacetime may provide a way out of the 
ultraviolet problem with the reduction of spacetime dimensions. It is inter- 
esting to note that in addition to these different approaches to quantum 
gravity which indicate the emergence of a ground-scale spacetime with frac- 
tal properties, there also exists a recent work 12 which shows that spacetime 
with quantum group symmetry in general has scale-dependent dimensions, 
and it allows one to establish a link between non-commutative geometries 
and quantum group and seemingly unrelated approaches of quantum grav- 
ity like causal dynamical triangulations and quantum Einstein gravity. All 
these different theories have a common property, namely, dimensions of 
spacetime depend on scale, and they could have fractal properties at small 
scales. At a more fundamental level, Palmer 13 introduced the postulate of 
an invariant set which assumes the existence of a state space embedding 
an invariant fractal subset which represents physical reality. It provides a 
possible link between fractal geometry and quantum theory and it may help 
to resolve some basic problems in quantum mechanics. 

One useful mathematical tool in dealing with fractals and fractal-related 
phenomena is fractional calculus 14-18 . An important development that 
stimulates the application of fractional calculus to physical sciences is the 
realization of the close connection between fractional calculus and fractal 
geometries 19-22 . Such a relation makes fractional differential equations a 
natural tool for describing transport processes with fractal properties, such 
as anomalous diffusion, non-Debye relaxation processes, and other frac- 
tal phenomena 16,17 ' 23,24 . Although the application of fractional differential 
equations to model various transport phenomena in complex heterogeneous 
media becomes more and more popular during the past two decades, their 
use in quantum theory is still quite limited. For many years the trajectories 
of the well-known fractal process, the Brownian motion, have been exten- 
sively used in Feynman path integral approach to quantum mechanics 25 . 
The paths of the Brownian motion are known to be fractal curves with 
Hausdorff dimension two 26 . Brownian motion also plays an important role 
in stochastic mechanics, which provides an alternative formulation to quan- 
tum mechanics 27 . In the early application of fractal geometry in quantum 
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field theory, the main focus was on the quantum field models in fractal sets, 
and quantum field theory of spin systems such as Ising spin model 28 ' 29 . 
Fractal geometry was later applied to Wilson loops in lattice gauge theory, 
and in the modeling of random surfaces in quantum gravity 30 . 

The possible existence of fractal structure of quantum spacctime 
prompted one to employ fractional calculus and fractional differential equa- 
tions 1 in quantum theory. Just like fractional diffusion equation has re- 
placed diffusion equation to describe anomalous transport systems in frac- 
tal media, for quantum theories in fractal spacctime it is natural to con- 
sider quantum mechanics and quantum field theory which satisfy fractional 
generalizations of Schrodingcr equation, Klein-Gordon equation and Dirac 
equation. There are some studies on the formulation of quantum mechanics 
based on various types of fractional Schrodingcr equations and their possi- 
ble applications 31-38 . By noting that fractional diffusion equations are just 
the Euclidean Schrodingcr equations, results obtained for the former can 
be directly applied to the latter. In particular, fractional Levy path integral 
approach of fractional quantum mechanics has been considered 31,32 ' 39,40 . 
It will be interesting to see whether it is possible to use the path integral 
representation of fractional Brownian motion 41,42 and fractional oscillator 
processes in the formulation of fractional quantum theory. It has been 
shown that fractional Brownian motion does appear naturally in the frac- 
tional gluon propagator with the temporal gauge condition 44 . Incorpora- 
tion of supcrsymmetry into fractional quantum mechanics has also been 
considered 45 . 

Fractional Klein-Gordon equation and fractional Dirac equation have 
been studied by several authors during the past decades 46-48 . Fractional 
power of D' Alembertian operator was used in the nonlocal kinetic terms La- 
grangian field theory in the (2+l)-dimensional bosonization. Such fractional 
power operators also appeared in the effective field theory which has some 
degrees of freedom integrated out from the underlying local theory 49-51 . 
Canonical quantization of fractional Klein-Gordon massless and massive 
fields has been studied by some authors 52 ' 53 . Quantization of fractional 
Klein-Gordon field and fractional gauge field based on Nelson's stochas- 
tic mechanics and Parisi-Wu stochastic quantization procedure at zero and 
positive temperature have been considered 44 ' 54 . There also exists work on 
axiomatic approach to fractional Klein-Gordon field, where properties of the 
n-point Schwinger or Euclidean Green functions and their analytic contin- 
uation to the corresponding n-point Wightman functions were given 55,56 . 
Fractal fermion propagator has been obtained as a consequence of the QED 
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radiative corrections, which requires the incorporation into the propagator 
with a fractional exponent connected with the fine structure constant 57 . 
More recent work is related to the Casimir effect for the massless and mas- 
sive fractional fields at zero and positive temperature; and work on Casimir 
effect due to fractional Klein-Gordon field subject to fractional Neumann 
boundary conditions which interpolate between the usual Dirichlet and 
Neumann conditions have been carried out 58 . 

In this paper, we shall mainly concern with the Casimir effect associated 
with fractional Klein-Gordon fields. In the next section a brief review on the 
results obtained for the Casimir energy and Casimir force due to a fractional 
Klein-Gordon field will be given. We shall also discuss the implication of 
fractional Neumann boundary conditions on the Casimir force. Calculations 
of Casimir force related to piston geometry for the massive fractional Klein- 
Gordon field will be given in detail in Section 5 . In the concluding section 
some possible directions for further work will be discussed. 

2. Casimir Effect Associated with Fractional Klein-Gordon 
Field 

Currently, the biggest challenge in cosmology is the explanation of the ac- 
celerated expansion of the universe. This is closely related to the problem 
of dark energy which has been proposed to explain the cause of this accel- 
erating expansion. Various candidates of dark energy have been considered, 
among them are positive cosmological constant, quintessence, Casimir en- 
ergy in extra spacetime dimensions, etc 59 ' 60 . As a possible candidate for 
dark energy, Casimir energy 61-64 has attracted considerable renewed in- 
terest. On the other hand, Casimir effect has become significant in nan- 
otechnology. At the nano-scales, the Casimir force cannot be neglected, 
and therefore scientists working in nanotechnology begin to show interest 
in this subject. Improvement in the ability to measure small forces near 
surfaces in nano-devices has allowed the possibility to manipulate Casimir 
force to drive a nano-device. Experimental advances together with improved 
computations and simulations of Casimir forces make it possible to tailor a 
nano-device to eliminate unwanted properties caused by the Casimir force 
such as stiction or adhesion, and hence facilitate the design of better and 
more effective nano-devices. Thus, from both large and small scales Casimir 
effect seems to play an important role, as it is expected from vacuum energy 
which is all prevalent in the universe 65 . 

In this section we summarize results on some physical quantities related 
to Casimir effect due to fractional Klein-Gordon fields. Denote by A = 
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Q2 u Q2 

—rr + — — 7? the (D + l)-dimensional Euclidean Laplacian operator. The 
ot t-r 1 oxj, 
j 3 

Riesz fractional derivative (— A) Q is defined as: 14,18 

(-Arf(t,x)=F- 1 ((k 2 + k 2 ) 2a f(k a ,k)) , 

where /(fco,fc) = F(f(t,x)) is the Fourier transform of f(t,x). Let 4>(t,x), 
t £ K, x £ M. D , be the real scalar Klein-Gordon field. Consider the La- 
grangian density 

L = --4>(t,x)k{d)<j,{t,x), 

where A(d) is a pseudo-differential operator. The resulting Lagrangian 
equation is then given by A(d)(j)(t, x) = 0. For the scalar massive frac- 
tional Klein-Gordon field, one can have the following possible choices of 
fractional Klein-Gordon operator: 

A(d) =(-A + m 2 y (I) 
A(d) =(-A) a + m 2a (II) 
A(d) = [{-A) a +m 2a Y (III) 

Here we remark that I seems to be the simplest choice, and II is a special 
case of III when 7 = 1. Simple results on Casimir effect associated with 
fractional Klein-Gordon field of type I have been obtained 54 . Note that 
the zero temperature fractional Klein-Gordon scalar massive field with Eu- 
clidean propagator is given by S(k) = (k 2 + m 2 ) -7 ; and the Euclidean 
propagator for the positive temperature case takes the following form: 
S{k n ) = (fc 2 + m 2 )' 7 , with k 2 = uj 2 + k 2 , u n = 2nir/P, and 1//3 = T 
is the temperature. The free energy corresponds to the type I fractional 
Klein-Gordon field equals to that for ordinary Klein-Gordon field multiply 
by a factor 7. When D = 3, the high temperature limit of the free energy 



is 54 



^ T 4 + l m 2 T 2_^ m 3 r . 

90 24 12tt 



In general, it is very difficult to compute the free energy for the type III 
fractional Klein-Gordon field. The only case it has been computed is when 
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Here Cr( s ) = n S ^ s the Riemann zeta function and u) a = 1 if and 

n=l 

only if l/(2a) is an integer. Otherwise ui a = 0. The dependence of the 
dimensionless free energy J 7 /m on /3m and a are shown in Fig. 1 and Fig. 
2. 




Fig. 1. T/m as a function of 0m when a = 0.1,0.3,0.5,0.7,0.9 and 7 = 1 



Casimir energy due to massless field between two parallel plates subject 
to fractional Neumann boundary conditions has been computed based on 
the zeta function regularization method 85 . Let us recall briefly the basic 
steps involved in this regularization procedure which can be summarized 
in three steps. In the case of scalar massless fractional Klein-Gordon field, 
they are (i) Evaluation of eigenvalues A of the fractional D'Alembertian 
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Fig. 2. F/m as a function of /3m and a. Here 7 = 1. 

operator (— A) Q with certain boundary conditions, thus leading to spec- 
tral zeta function C(-A)°( s ) = ^""^ A~ s ; (ii) Analytic continuation of the 

A 

zeta function to a meromorphic function of the whole complex plane; 
and (iii) Calculation of dct(— A) a in terms of C(-A)°(s) using the relation 

det(-A) Q =exp(-C ( '_ A)Q (0)). 

Fractional Neumann boundary conditions were introduced to study the 
Casimir effect due to fractional Klein-Gordon field. On two plates separated 
by a distance a, the fractional Neumann boundary conditions are: 



— <f>{t,x,x D ) 



dx D 



— <f)(t,X,X Dj 



dx D 



= 0, 

x o=0 

= 0, 



with x G R-° 1 , X£> E K and /i, v G [0, 1], and the fractional derivative is 
defined in terms of its Fourier transform: 
d»f(x) 1 



dk{ikYe™ x f{-k), 



where (±ik) a = |fc| a e ±iQ7r / 2 sign(fc). 

Results on free energy density and pressure due to fractional massless 
field confined to parallel plates subject to the following fractional Neumann 
boundary conditions have been obtained 58 . These include (a) General frac- 
tional Neumann boundary conditions with p ^ f and v) ^ (0, 1) or 
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(1,0); (b) \x = v 7^ 0,1; (c) Dirichlet boundary conditions [i = v = 0; 
(d) Neumann boundary conditions ji = v = 1; and (e) Boyer boundary 
conditions (/i, v) = (0, 1) or (1, 0). Low and high temperature limits of the 
Casimir energy and pressure have been calculated for the above cases. It is 
interesting to note that there exists a transition value in the difference in 
the orders of the fractional Neumann conditions for the two plates, that is 
|/i — v\, for which the Casimir force changes from attractive to repulsive, or 
vice versa. 



3. Topological Symmetry Breaking of Self-interacting 
Fractional Klein-Gordon Field 

Quartic self-interacting fractional Klein-Gordon scalar field in a toroidal 
spacetime and the problem of topological mass generation and symmetry 
breaking has been investigated 67 . The Lagrangian is 

L = -^<f>(t,x)A(d)cf>(t,x) - ^(t,x), 

where A(<9) = (—A + ?ti 2 ) 7 is the type I fractional Klein-Gordon operator. 
Based on the Epstein zeta function regularization method, the one-loop ef- 
fective potential for both the massive and massless </> 4 theory in the toroidal 
spacetime T p x R 9 is derived in terms of power series of \(p 2 , where (p is the 
constant classical background field. Renormalization of these quantities is 
carried out, and the results for the renormalized mass are obtained explic- 
itly. In the massive case, the renormalized mass m ren is related to the bare 
mass m by 



,2 7 =TT1 2 7 | 

(27r)i+T( 7 ) 



m ten = m 7 + " — — X 



P 



d-2j_ 



kgZP\{0} \i=l / 



Here d = D + 1 = p + q, Li, . . . ,L p are the compactification lengths of 
the torus T p . In the massless m = case, the theory is nonrenormalizable 
when a = d/2. When a ^ d/2, the renormalized topologically generated 
mass m ron is 

• If p = 0, then rr&n = 0; 



If p > 1 , and 
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Here Ze, p (s; 01, . . . , a p ) is the zeta function 



2tt 




2tt\ 



Ze,p(s; 01, . . . , a p ) 



From these, conditions for symmetry breaking are derived analytically, 
which show that there is no symmetry breaking in the massive case, how- 
ever, for the massless case, symmetry breaking occurs provided < d— 2j < 
p for a fixed compactified dimension p < 9. In the case p > 10, symmetry 
breaking can only appear for d — 27 in a proper subset of (0,p\. 

In Fig. 3, Fig. 4, Fig. 5, Fig. 6 and Fig. 7, the regions where the renor- 
malized topologically generated mass m 2 ^ n are positive and negative respec- 
tively are shown. 



4. Casimir Piston Associated with Massless Fractional 
Klein-Gordon Field 

In the calculations of Casimir force based on the conventional methods for 
a confined geometry such as a rectangular cavity, one usually ignores the 
nontrivial contribution of vacuum energy outside the cavity. A finite re- 
sult is obtained by discarding the surface divergent terms which depend 
on the dimensions and geometry under the pretext of some regularization 
methods 61-64 . However, there are studies which show that renormalization 
of physical parameters of the theory fails to remove such surface diver- 
gence 68 ' 69 . 

In 2004, Cavalcanti 70 introduced a new geometric setup, namely a two- 
dimensional rectangular piston and showed that for such a geometrical 
configuration it is possible to obtain an unambiguous finite Casimir force 
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Fig. 3. The regions where m r J„ > and m T J n < for p = 2 and L1L2 = 1. Here 




0.2 0.4 0.6 0.8 1 1.2 1.4 



s 

Fig. 4. The regions where m^2n > and m^Jn < for p = 3, L1L2L3 = 1, L\ : L2 : 
£3 = k : 1 : 1. Here 8=5—7. 

despite the presence of the divergence in the Casimir energy. He demon- 
strated that in the case of a massless scalar field with Dirichlct boundary 
conditions, a finite attractive Casimir force without the surface divergence 
is obtained due to the cancelation of the divergent part of the Casimir 
force acting in the two regions separated by the piston. This result leads 
to a surge in the studying of the Casimir force on pistons with various 
geometric configurations and boundary conditions 71-87 . Casimir force on 
a rectangular piston in the spacetime with extra compactified dimensions 
has also been studied 88-93 . It has been shown that the Casimir force be- 
tween two (nonmagnetic) dielectric bodies which are related by reflection 
is always attractive 94,95 . The attractive nature of the Casimir force may 
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-2-1 1 2 

Fig. 5. The region where rn^Jn > and m^2 n < for p = 3, L1L2L3 = 1, L\ : L2 
Li = 1 : fc 2 : k 3 and f - 7 = 0.1. 
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logk 



Fig. 6. The regions where m^2n > and m^en < for p = 3, L1L2L3 = 1, £1 : £2 : 



L3 = 1 : k,2 : fcsand 



"7 : 

2r 
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_ 2 i , , , 

—2 — 1 1 2 

l°gk 2 

Fig. 7. The regions where m^n > and < for p = 3, L1L2L3 = \, L\ : L2 '■ 

L 3 = 1 : fc 2 : ^3 and f - 7 = 0.6. 
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result in undesirable effects stiction in nano devices. As a result, it is desir- 
able to search for circumstances under which the Casimir force can become 
less attractive, or even repulsive. Results obtained by Barton 77 showed that 
for a thin piston with weakly reflecting dielectrics, the nature of Casimir 
force is separation-dependent. The force is attractive at small separations, 
but it becomes repulsive when the separation increases. For Casimir ef- 
fect on an n-dimensional (n = 1, 2 or 3) rectangular piston with mixed 
boundary conditions, namely with one surface having Dirichlct condition, 
and the other with Neumann condition, the resulting force is repulsive. 78 
Other possible ways of obtaining a repulsive Casimir force have also been 
discussed 79 ' 83 ' 85 ' 86 ' 92 ' 96 . In particular, the Casimir force of fractional Klein- 
Gordon field acting on a piston with fractional boundary conditions can be 
repulsive 85,96 . 

Finite temperature Casimir force acting on a (D+ l)-dimcnsional piston 
in a (D + l)-dimcnsional semi- infinite Dirichlct cylinder due to a type I 
fractional masslcss Klein-Gordon field at positive temperature has been 
calculated 85 ' 96 . For such a system, when the fractional Neumann condition 
of order [i is imposed on the piston, the Casimir force acting on the piston 
is given by 



where a is the distance between the piston and the opposite wall; {rrij} is 
the set of eigenvalues of the Laplacian operator with Dirichlct boundary 
conditions on the cross section of the cylinder. 

This result has to be divided by two if fj, = (Dirichlet) or \i = 1 
(Neumann). From this, we find that for fi > 1/2, the Casimir force acting 
on the piston is always repulsive, and the magnitude of the force decreases 
with a. On the other hand, when /i < 1/2, the Casimir force can change 
from repulsive to attractive with increasing a. For certain values of /x less 
than but close to 1/2, the Casimir force can be attractive or repulsive 
depending on the aspect ratio of the cavity and the temperature. Increasing 
temperature can change the Casimir force from attractive to repulsive. In 
the case D = 3 and the cross section of the cylinder is the rectangle box 
[0,^2] x [0, L3], the dependence of the Casimir force on various quantities 
are shown graphically in Fig. 8, Fig. 9, Fig. 10 and Fig. 11. 
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Fig. 8. The Casimir force F as a function of < 
and T = IK, 2K and 5K respectively. 



when d = 3, L 2 = L 3 = 0.01m, /j, = 0.47, 




Fig. 9. The Casimir force F as a function of /i when d = 3, L2 = ^3 = 0.01m for 
various values of (a, T) . 



5. Casimir Piston Associated with Massive Fractional 
Klein-Gordon Field 

In this section, we present some new results on the Casimir effect of a 
massive fractional Klein-Gordon field in the piston scenario. Consider the 
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0.5 1 1.5 2 2.5 3 

T(K) 

Fig. 10. The Casimir force F as a function of T when d = 3, L2 = L3 = 0.01m, fi = 0.47 
and a = 0.8mm, 0.9mm, 1mm and 1.2mm respectively. 




-400 1 1 1 1 ' 1 ^ 
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T(K) 



Fig. 11. The Casimir force F as a function of T when d = 3, L2 = Lz = 0.01m, fi = 0.48 
and a = 0.8mm, 0.9mm, 1mm and 1.2mm respectively. 



type III fractional Klein Gordon field (f>(t, x) satisfying the equation: 



[(-Ar + m 2 f #,x)=0. 
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We are interested in the finite temperature Casimir effect of the fractional 
Klein-Gordon held <fi(t, x) on a piston inside a semi-inhnite rectangular 
cylinder [0,oo] x [0,^2] x . . . x [0,Ld], with Dirichlct boundary conditions 
on the walls of the cylinder and fractional Neumann condition of order [i on 
the piston. The rectangular piston is located at x\ = a. To regularize the 
Casimir energy, hrst assume that the system is enclosed in the rectangular 
box [0, Li] x ... x [0, Ld\- Take another piston system where the piston is 
located at x\ = ii/77, for a constant 77 > 1. The regularized Casimir energy 
is defined as the limit (see Fig. 12) 



system A 



A 


A A 






i-jj/ 













Fig. 12. Regularizing of Casimir energy. 



where E^ as and E^ as are respectively the Casimir energies of system with 
the piston at x± = a and at x\ = L\jr\ respectively. 



-^Cas = ^Cas( a )+^Cas(-kl- a )> -^Cas = E Ca.s 



Ll 



vJ Cas V 1 v 



Using zeta regularization, the finite temperature Casimir energy inside the 
box [0, a] x [0, L 2 ] x . . . x [0, Lb] is given by 



£c°as(«) = -f (&(0| «) + log[^ 2 ]CT(0; a)) , 
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where v is a normalization constant, and (r(s;a) is the zeta function 
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When /i = or [i = 1, the summation over n starts from n = 1. Using 
standard techniques, 
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where 

oo 

k^eN 15 - 1 ;=o 

Therefore, 

/■ i \ vi \ , 8q ^ (-iy 2 , r(7s + j) ^ \ - 

Ct(s; a) =aY(s) H =— — -> — - — m J _, . . > > > 

V V' ^ J= J y w J; ' n=lk i £N°- 1 1=0 

where i^(s) is independent of a. One immediately finds that 

C T (0;a) = aY 1 (0), 

and 

V j=0 V J ; n=l k^eN " 1 Z=0 

(2naa; kX)i ). 

From this, we find that the finite temperature Casimir energy in the box 
[0, a] x [0, L 2 ] x ... x [0, L D ] is given by 

V j=0 V 7 'n^lkxgNfl-itO 

X K aj-i ( 2nauJ k ± ,l) , 

where the term i?i is independent of a. This term will vanish when we 
take the difference of the energies in system A and system B. In the limit 
L\ — > oo, we find that the finite temperature Casimir energy of the piston 
system is 

*-~ ^E(-)'F7^T)E E E' e """ 

V J=0 V J ; n=lk i £N D - 1 i=0 

/ na \ QJ ~5 
x Jf ._i (2naw ki , ; )- 

For /j = or \i = 1, this result should be divided by two. Eq. (1) is the 
high temperature expansion of the Casimir energy. It shows that the high 
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temperature leading term is given by the sum of the terms with I = 0, i.e. 
when aT ^> 1, 

V ]=0 V J > n=lk ± 6N D - 1 

/ na \ QJ "5 



The zero temperature Casimir energy is obtained by taking the limit T — > 
of (1), which gives 



ttI ^ ( ' T(aj + 1) ^ ^ 



^Cas 



<\l- 



oo 



OO OrvT 00 / \ Oil — 1 

j=0 v J ' „=lk 1 £N 1 '- 1 v x 7 

x (2naw ki ) . 

The Casimir force acting on the piston is given by 

,11 _ d4, s 



^Cas 



da 

8 Q7 T V ■ m J v ^ \ - 

v j=o „=i ki eN" 



i^cj+i (2naw k _ L , z ) 



v j = l v J; n = l lc^eN 13 " 1 Z=0 



;/(/ 



\w kj 



# ■_ i (2naa; kxi • 
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In the limit T — >• 0, wc find that the zero temperature Casimir force acting 
on the piston is given by 



„,T=o_ 2a 7>r ^ ^ 

Cas ~ TT ^ [ ' T(aj + 1) L> ^ 



OO oo 



S^(-l)i(2ao-l) = 



7r * — ' ' r (aj + 1) ^— ' ^ 

y J 1 n=lk, PN"" 1 



j=o v J y n=i kj_eF 

aj— 1 

ifcy-i (2naw ki ) 



- 



no 

For a type III massive Klein-Gordon field, we can only expand the Casimir 
energy and Casimir force as a series in m 2a . It is not easy to determine the 
sign of the Casimir force from these series expressions. 



6. Conclusion 

Justifications on the need to study fractional Klein-Gordon fields are given. 
A brief summary of the results on fractional Klein-Gordon fields mainly by 
the author and coworkers are given. We also present new results on Casimir 
piston associated with fractional Klein-Gordon massive field. 

Finally, we would like to suggest some other possible directions for fur- 
ther work. So far we have considered fractional Klein-Gordon field with 
unique fractional exponent. Even though type III fractional Klein-Gordon 
field parametrized by two indices provides more flexibility for modeling, it 
is still inadequate to represent quantum field in spacetime with dimensions 
that vary with spacetime scales. In view of the possible variable spacetime 
dimensions, one may have to deal with fractional quantum fields of variable 
fractional order, with replaced by the function . Fractional stochastic pro- 
cesses of variable order such as multifractional Brownian motion, 97,98 mul- 
tifractional Levy motion, 99 ' 100 multifractional Riesz-Bessel process 101 and 
multifractional Ornstein-Uhlenbeck process 102 have been studied. Variable 
order one-dimensional Klein-Gordon massive field is given by multifrac- 
tional Ornstein-Uhlenbeck process, which can be easily extended to higher 
dimensions to yield multifractional Klein-Gordon field. On the other hand, 
fractional Klein-Gordon field can also be obtained as a special case of mul- 
tifractional Riesz-Bessel field. In addition, it may be necessary to impose 
fractional Neumann boundary conditions of variable fractional order when 
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dealing with multifractional Klein-Gordon field. It should be emphasized 
that calculations of physical quantities of Casimir effect due to fractional 
Klein-Gordon field of variable order arc expected to be mathematically dif- 
ficult, and new techniques will be required. 

Another direction of generalization is to consider the above results 
in non-flat space, which requires fractional itegro-differential operators in 
curved spacetime. One particular case of recent interest will be to study 
fractional Klein-Gordon field in de Sitter spacetime. Stochastic quantiza- 
tion for ordinary Klein-Gordon field in dc Sitter spacetime has been carried 
out. 103 It will be interesting to extend it to the fractional case. Path inte- 
gral approach to fractional Klein-Gordon field is another area that has so 
far not been seriously treated. 
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